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We study two different types of simplified models for Kelvin wave turbulence on quantized vortex 
lines in superfluids near zero temperature. Our first model is obtained from a truncated expansion 
of the Local Induction Approximation (Truncated-LIA) and it is shown to possess the same scalings 
and the essential behaviour as the full Biot-Savart model, being much simpler than the latter 
and, therefore, more amenable to theoretical and numerical investigations. The Truncated-LIA 
model supports six-wave interactions and dual cascades, which are clearly demonstrated via the 
direct numerical simulation of this model in the present paper. In particular, our simulations 
confirm presence of the weak turbulence regime and the theoretically predicted spectra for the 
direct energy cascade and the inverse wave action cascade. The second type of model we study, the 
Differential Approximation Model (DAM) , takes a further drastic simplification by assuming locality 
of interactions in fc-space via a differential closure that preserves the main scalings of the Kelvin 
wave dynamics. DAMs are even more amenable to study and they form a useful tool by providing 
simple analytical solutions in the cases when extra physical effects are present, e.g. forcing by 
reconnections, friction dissipation and phonon radiation. We study these models numerically and 
test their theoretical predictions, in particular the formation of the stationary spectra, and the 
closeness of the numerics for the higher-order DAM to the analytical predictions for the lower-order 
DAM . 

PACS numbers: 67.25.dk, 67.85.De, 47.37.+q 

INTRODUCTION 

It is well known that a classical vortex filament can support linear waves. These were predicted by Kelvin more than 
one century ago and experimentally observed about 50 years ago in superfluid He. At very low temperature, where 
the friction induced by normal fluid component can be neglected, Kelvin waves can be dissipated only at very high 
frequencies by phonon emission pQ . Therefore at lower frequency, energy is transferred among different wavenumbers 
by nonlinear coupling. This is the mechanism at the basis of the Kelvin wave cascade which sustains superfluid 
turbulence [2J [3] . 

In recent years, the single vortex Kelvin wave cascade has attracted much theoretical |H|S], numerical [S1|71|S] and 
experimental 10J attention. Even within the classical one-dimensional vortex model, different degrees of simplification 
are possible. For small amplitudes, the vortex configuration can be described by a two component vector field, made of 
the coordinates of the vortex line in the plane transverse to the direction of the unperturbed filament. These depend on 
the single coordinate that runs along the filament. As was shown in [2J, this system of equations admits a Hamiltonian 
formulation, dubbed the two-dimensional Biot-Savart formulation (2D-BS), see (2) below. Another, more drastic, 
simplification is obtained by considering local interactions only. This leads to the local induction approximation 
(LIA) which was originally derived starting from the full 3D-BS [TT]. The main limitation of LIA is that it generates 
an integrable system with infinite conserved quantities, as it is equivalent to the nonlinear Schrodinger equation |12j . 
and therefore, the resonant wave interactions are absent (at all orders) and one cannot reproduce the phenomenology 
of the full system. For this reason LIA, despite its simplicity, is of little help for the study of weak Kelvin wave 
turbulence. On the other hand, LIA contains solutions leading to self-crossings (numerical [5] and analytical [5]) and, 
therefore, it can qualitatively describe vortex line reconnections in strong 3D-BS turbulence ("vortex tangle"). 

In this paper we consider simple models for a vortex filament that is able to sustain a turbulent energy cascade. 
The first model is obtained in the limit of small amplitudes by a Taylor expansion of the 2D-LIA. The truncation 
breaks the integrability of the Hamiltonian and therefore generates a dynamical system with two inviscid invariants 
(energy and wave action). For this class of systems, whose prototype is the two-dimensional Navier-Stokes turbulence 
|13j . we expect a dual cascade phenomenology in which one quantity flows to small scales generating a direct cascade 
while the other goes to larger scales producing an inverse cascade. The possibility of a dual cascade scenario for 
Kelvin waves turbulence has been recently suggested [5j[14] but never observed, and in this paper we present the first 
numerical evidence for the inverse cascade. 
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The second class of simplified models, Differential Approximation Models (DAMs), use a closure in which the 
multi-dimensional fc-space integral in the wave interaction term (collision integral in the wave kinetic equation) is 
replaced by a nonlinear differential term that preserve the main properties and scalings of the Kelvin wave dynamics 
such as the energy and wave action conservations, scaling of the characteristic evolution time with respect to the 
wave intensity and the wavenumber k. DAMs have proved to be a very useful tool in the analysis of fluid dynamical 
and wave turbulence in the past [5J [T5J [TBJ [T71 [TH1 [THl [23 121], and here we study them in the context of the Kelvin 
wave turbulence. DAMs are particularly useful when one would like to understand the temporal evolution of the 
spectrum, when the physical forcing and dissipation need to be included, or when the Kelvin wave system is subject 
to more involved boundary conditions leading to simultaneous presence of two cascades in the same range of scales, 
or a thermalization (bottleneck) spectrum accumulation near a flux-reflecting boundary . In the second part of this 
paper, we will present numerical studies of DAMs in presence of some of these physical factors and we will test some 
previously obtained analytical predictions. 



BSE 



At a macroscopic level, the superfluid vortex filament is a classical object whose dynamics is often described by the 
Biot-Savart equation (BSE) 

k f ds x (r — s) 

r = 4nJ |r-s| 3 ' W 

which describes the self-interaction of vortex elements. The quantum nature of the phenomenon is encoded in the 
discreteness of circulation n~h/m [2"2"] . 

The BSE dynamics of the vortex filament admits a Hamiltonian formulation under a simple geometrical constraint: 
the position r of the vortex is represented in a two-dimensional parametric form as r = (x{z),y{z), z), where 2 is a 
given axis. From a geometrical point of view, this corresponds to small perturbations with respect to the straight line 
configuration, i.e. the vortex cannot form folds in order to preserve the single- valuedness of the x and y functions. In 
terms of the complex canonical coordinate w(z,t) — x{z,t) + iy(z,t), the BSE can be written in a Hamiltonian form 
iw = SH[w]/6w* with [2, 

H[w] = ± f d z ld z 2 , l+^t^M , (2) 
4W ^/(z 1 ~z 2 )^ + \w(z 1 )~w(z 2 )\^ 

where we have used the notation w'{z) = dw/dz. The geometrical constraint of a small amplitude perturbation can 
be expressed in terms of a parameter e(zi, z 2 ) — \w(z\) — w(z 2 )\/\zi — z 2 \ <C 1. 

An enormous simplification, both for theoretical and numerical purposes, is obtained by means of the so called local 
induction approximation (LIA) This approximation is justified by the observation that ([TJ is divergent as s — > r 
and is obtained by introducing a cutoff at a < |r — s\ in the integral in (JT|) which represents the vortex filament radius. 

When applied to Hamiltonian ([2]), the LIA procedure gives [2] 

H[w] = 2^ In f ~J J dz^l + \w'(z)\ 2 = 2f3L[w], (3) 

where £ is a length of the order of the curvature radius (or inter- vortex distance when the considered vortex filament 
is a part a vortex tangle), [3 — (n / An) ln(£ / a) . Here, it was taken into account that because a is much smaller 
than any other characteristic size in the system, f} will be about the same whatever characteristic scale I we take 
in its definition. We remark that in the LIA approximation the Hamiltonian is proportional to the vortex length 
L[u>] = J dzyjl + |w'(z)| 2 which is therefore a conserved quantity. The equation of motion from |3| is (we set /? = 1/2 
without loss of generality, i.e. we rescale time as 2j3t — > t) 
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2 Vv^ 



(4) 



As a consequence of the invariance under phase transformations, equation Q also conserves the total wave action 
(also called the kelvon number [1]) 

N[w] = [ dz\w\ 2 . (5) 
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In addition to these two conserved quantities, the 2D-LIA model possesses an infinite set of invariants and is 
integrable, as it is the LIA of BSE (which can be transformed into the nonlinear Schrodinger equation by the Hasimoto 
transformation, see Appendix B). Due to the integrability, in weak Kelvin wave turbulence [30] the energy and the 
wave action cannot cascade within the LIA model, but this can be fixed by a simple truncation as we show in the 
next section. 



TRUNCATED LIA 

Integrability is broken if one considers a truncated expansion of the Hamiltonian ^ in power of wave amplitude 
w'(z). Taking into account the lower order terms only, one obtains: 

H exp [w] = H + H 1 +H 2 = 

= J dz(i+ l -\ w r-\w\^. (6) 

Neglecting the constant term, the Hamiltonian can be written in Fourier space as 

H eX p = J UJ k \w k \ 2 dk+ J dk 12 34W 1 234$llw* 1 W2W 3 W 4 , (7) 

with ui — k 2 /2, W1234 = — |fcifc2&3fc4 an d we used the standard notation <5|f = 5(ki + k 2 — k 3 — foj) and dfci234 = 
dkidk 2 dk 3 dk4 : . 

In Wave Turbulence, the near-identity transformation allows one to eliminate "unnecessary" lower orders of nonlin- 
earity in the system if the corresponding order of the wave interaction process is nil |23j . For example, if there are no 
three-wave resonances, then one can eliminate the cubic Hamiltonian. (The quadratic Hamiltonian corresponding to 
the linear dynamics, of course, stays). This process can be repeated recursively, in a way similar to the KAM theory, 
until the lowest order of the non-trivial resonances is reached. If no such resonances appear in any order, one has an 
integrable system. 

In our case, there are no four-wave resonances (there are no non-trivial solution for the resonance conditions for 
lu ~ k x if x > 1 in one dimension). However, there are nontrivial solutions of the six- wave resonant conditions. 
Thus, one can use the near-identity transformation to convert system ^ into the one with the lowest order nonlinear 
interaction to be of degree six, (there are no five-wave resonances, the interaction coefficients in the quintic Hamiltonian 
are identically equal to zero after applying the canonical transformation.) 

A trick for finding a shortcut derivation of such a transformation is described in |23j . It relies on the fact that the 
time evolution operator is a canonical transformation. Taking the Taylor expansion of w{k, t) around w{k, 0) = c(fc, 0) 
we get a desired transformation, that is by its derivation, canonical. The coefficients of each term can be calculated 
from an auxiliary Hamiltonian H aux , 
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Vl23#3 2 (C1C2C3 + C.C.) dfci23 

/ U123S 123 (c lC2 c 3 + c.c.) dk x 



W^1234^34 c l c 2 c 3C4rffc 



123 L 



J -^1234^234 {ci c *2 c *3 c X + c - c -) *1234 



+ / Yl234^ 1234 (C1C2C3C4 + C.C.) dfc 



1234 



J ■Z'l2345^ 12345 (C1C2C 3 C 4 C 5 + C.C.) dfcl2345 
J ^12345^2345 ( C 1 C 2 C 3 C 4 C 5 + c - c ^12345 
J -Bi2345^345 ( c lC2C3C4 c 5 + c - c *12345 
J C*i23456^456 c l c 2 c 3 c 4 c 5 c 6 'fcl23456- 



(8) 



The auxiliary Hamiltonian H aux represents a generic Hamiltonian for the canonical variable Cfc, thus, defining c& 
in the canonical transformation will set the interaction coefficients of H aux . Here all interaction coefficients, (terms 
denoted with tildes) present in the auxiliary Hamiltonian are arbitrary. A similar procedure was done in Appendix 
A3 of [23 to eliminate the cubic Hamiltonian in cases when the three-wave interaction is nil, and here we apply a 
similar approach to eliminate the quadric Hamiltonian. The transformation is represented as 

The transformation is canonical for all t, so for simplicity we set t = 1. The coefficients of Q can be calculated from 
the following formulae, 

dc(k,t) \ .SHaux 

dt ) t=0 - 1 5c* ' 
d 2 c(k,t)\ .d6H aux 

-BP-)^ = ^dt^*-- (10) 
Due to the original Hamiltonian (oj H exp having U(l) gauge symmetry, we have no cubic order Hamiltonian terms, 



this greatly simplifies the canonica 



transformation Q9J) , because the absence of any non-zero three- wave interaction 
coefficients in H exp automatically fixes the arbitrary cubic (and quintic) interaction coefficients within the auxiliary 
Hamiltonian H aux to zero. Thus, transformation (fOp reduces to 



Wfc = Ck ~ ^ J dk2SiWk234C2C3C^l - 3i /'dfc23456C , fe23456^456 c 2 c 3 C 'l c 5C6 

+ \ J &234567(w^743HW4 fe 5^2 7 3 6 " 2Wfc247 WW*?^) C* 2 C* 3 C 4 C 5 C 6 . (11) 

To eliminate the nonresonant four- wave interactions in Hamiltonian Q, we substitute transformation (111 into 
Hamiltonian (7). This yields a new representation of Hamiltonian (7) in variable where the nonresonant terms 
(more specifically the four- wave interaction terms) will involve both W1234 and W\ 23 a- Arbitrariness of W1234 enables 
us to select this to eliminate the total four- wave interaction term, H 2 . In our case this selection is, 

^1234 = + Wl2 " • (12) 

U>i + uo 2 — 0J3 — OJ4 
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This choice is valid as the denominator will not vanish due to the nonresonance of four- wave interactions. Hamiltonian 
([7]) expressed in variable Ck, H c becomes 



H c = J io k c k c* k dk 

+ J (Cl23456 - i (uJl + U>2 + U 3 - UJi ~ UJ 5 - CJ 6 ) 

X C*i23456 ) ^456 c i c 2 c 3 c 4C5C6*i23456- (13) 



C123456 is the arbitrary six-wave interaction coefficient arising from the auxiliary Hamiltonian. This term does 
not contribute to the six-wave resonant dynamics as the factor in front will vanish upon the resonant manifold, 
that appears in the kinetic equation. C123456 is the six- wave interaction coefficient resulting from the canonical 



transformation which is defined later in equation (14) 



To deal with the arbitrary interaction coefficient C123456, one can decompose C123456 into its value taken on the 
six-wave resonant manifold plus the residue (i.e. C123456 = — ^fcifefcs^fcs&e + Ci23456)- We can then choose 
C123456 = — *Ci23456/( tJ i + w 2 + w 3 — w 4 — ^5 ~~ ^6); hence, allowing the arbitrary six-wave interaction coefficient in 
H aU x to directly cancel with the residual value of Ci23456- This enables us to write Hamiltonian H c as 



H c = J UJk\c k \ 2 dk + J dfci23456Ci23456<5456 c i c 2 c 3 c 4C5C6, 

The explicit form of the interaction coefficient C123456 can be expressed by 



C 



123456 



1 
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E 



E 



w j+k - rrJk 

(u>j + k - r +UJ r —Uij—UJ k ) 

Wi + j- pp ijW q+T - kkqr 
{uj q + r _ k +u> k -u> q -uj r ) ■ 



(14) 



Zakharov and Schulman discovered a parametrisation |24j for the six- wave resonant condition with ujk ~ fc 2 , 



h 

h 



P + R 
P + R 



1 1 

u v 
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- + - 

U V 



3v 
3v 
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u 

P+ — -2Ru, 
u 



(15) 
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P+R 



1 1 
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U V 

1 1 

u 

U V 



3v 
3v 



This parametrisation allows us to explicitly calculate C123456 upon the resonant manifold. This is important because 
the wave kinetics take place on this manifold, that corresponds to the delta functions of wavenumbers k and frequencies 
ui k within the kinetic equation. When this parametrisation is used with equation (14) and W1234 = —^1/02^3^4 we 



find that the resonant six-wave interaction coefficient simplifies to C123456 = — y§ k\ ki k% k^ k^ /c 6 . Note, that this is 
indeed the identical to the next term, H 3 in the LIA expansion with opposite sign. 

This six-order interaction coefficient is obtained from coupling of two fourth-order vertices of H 2 ■ It is not surprising 
that the resulting expression coincides, with the opposite sign, with the interaction coefficient of H3 in ([3]). Indeed, 
([3]) is an integrable model which implies that if we retained the next order too, i.e. H3, then the resulting six-wave 
process would be nil, and the leading order would be an eight-wave process in this case. In fact, in the coordinate 
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space the Hamiltonian is simply 



H exp [w] = H 



dz 



H 3 
1 



,/|2 



„'|6 



16' 



(16) 



Thus, the existence of the six- wave process is a consequence of the truncation ^ of the Hamiltonian. 

Hamiltonian Q (or equivalently (14) or (16l) constitutes the truncated-LIA model for Kelvin wave turbulence. It 
possesses the same scaling properties as the BSE system: it conserves the energy and the wave action, and gives rise 
to a dual-cascade six-wave system with an interaction coefficient with the same order of homogeneity as the one of 
the BSE. A slight further modification should be made in the time re-scaling factor as (3 — (n/Air), - i.e. by dropping 
the large log factor from the original definition. 

Physical insight into Kelvin wave turbulence is obtained from the wave turbulence (WT) approach which yields a 
kinetic equation which describes the dynamics of the wave action density n k — (|cfe| 2 ). 



The dynamical equation for the variable Ck can be derived from Hamiltonian ( 14) by the relation idck/dt = 5H c /5c* kl 
and is 



.dc k 
at 



= / dk 



2:iir,(> 



Cfe23456C2C3C4C5C6(54|| 



(17) 



Multiplying equation (17 1 by c* k , subtracting the complex conjugate and averaging we arrive at 



- ^ — 6Jm ( I ^23456^23456 A-23456^456 ) 7 



(18) 



Where Jfe23456^|| = (4 C 2 C 3 C 4 C 5C 6 )- 

Assuming a Gaussian wave field, one can take J/C23456 to the zeroth order J^23456' which is simplified via Gaussian 
statistics to a product of three pair correlators, 



J, 



(o) 

fc23456 



n 2 n 3 n 4 [$f (<^tf| 
SI (5151 + 5*51) 



5§ (5j5! 



(19) 



However, due to the symmetry of (7^23456 this makes the right hand side of kinetic equation (21 1 zero. To find a 



nontrivial answer we need to obtain a first order addition J^3456 ^° ^fc23456- 



To calculate j\^J 34:56 one takes the 
time derivative of J/t23456i using the equation of motion (17) of the canonical variable Ck, and insert the zeroth order 



approximation for the tenth correlation function (this is similar to equation ( 19 1, but a product of five pair correlators 
involving ten wavevectors) . ^23456 can then be written as 



7 (1) - Bp iAu)t 

ty fc23456 — oe 



-4 



£23456 

Aw : 



(20) 



where Aw = u>k + Wi + W3 — W4 — 
^4^23456 = 3C^ 23456 nfcn2n3n4n5n 6 



W 5 - 



loq and 
f ^ + ^ 

"2 "3 



J_ 

ri4 



The first term of (20 1 is a fast oscillating function, 
its contribution to the integral (|18| decreases with z and is negligible at z larger than 1 /w k , and as a result we will 



ignore the contribution arising from this term. The second term is substituted back into equation (18 1, the relation 
Jm(Aw) ~ — t:5(Aw) is applied because of integration around the pole, and the kinetic equation is derived, 



h k — I8n I rffc 2 3456 I Cfc 2 3456 1 2 ^456 ^( w 456) /fc23456, 



(21) 



where we have introduced /fc23456 = n k n-2,nzn 4 n^rie 
w 4 - w 5 - we)- 



nfc n-2 ri3 714 



and 5(w^) = 5(w k + w 2 + w 3 



A simple dimensional analysis of (21 1 gives 



n k ~ k 1 n k 



(22) 
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which is the same form obtained from the full BSE [3]. 

In wave turbulence theory, one is concerned with non-equilibrium steady state solutions of the kinetic equation 
(21 1. These solutions, that rely on a constant (non-zero) flux in some inertial range are known as Kolmogorov- 



Zakharov (KZ) solutions. In addition, the kinetic equation (211 contains additional solutions that correspond to the 
thcrmodynamical cquipartition of energy and wave action. These equilibrium solutions stem from the limiting cases 
of the more generalised Rayleigh- Jeans distribution 



T 



n k = 



(23) 



where T is the temperature of the system and \x is a chemical potential. 

To find the KZ solutions, one can apply a dimensional argument on both energy and wave action fluxes. The energy 
flux at wavenumber k is defined as III! = f dk'h^iu>k> which, using (22), becomes Ilj^ ~ fc 17 n|. By requiring the 



existence of a range of scales in which the energy flux is fc-independent leads to the spectrum 

fc- 17 / 5 , 



which, again, is the same form obtained from the full BSE 

A similar argument can be applied to the wave action (15 ) whose flux is Il fe 
independent flux of wave action requires a spectrum [M 



(AT) 



J dk'h k i 



(24) 



Therefore a scale 



k- 



(25) 



A word of caution is due about both spectra ( 24 I and ( 25 I because the dimensional analysis does not actually guarantee 



that they are true solutions of the kinetic equation. To check if these spectra are real solutions (and therefore physically 
relevant) one has to prove their locality i.e. convergence of the kinetic equation integral on these spectra. This has 
not been done yet, neither for the full BSE nor for the truncated-LIA model, and this is especially worrying since the 



spectrum ( 24 I has already been accepted by sizable part of the quantum turbulence community and has been used in 



further theoretical constructions. We announce that there is a work in progress to check locality of spectra (24 1 and 



(251 in both the BSE and truncated-LIA settings. However, as we will sec later, at least for the truncated-LIA these 



spectra are observed numerically, so we will tentatively assume that they are true and relevant solution. 



The two spectra ( 24 ) and ( 25 I occur in different scale ranges and the two cascades develop in opposite directions, 
as in the case of two-dimensional turbulence 13J. Among the two conserved quantities, the largest contribution to 
energy comes from smaller scales than those that contribute to wave action (because the former contains the field 
derivatives). Therefore, according to the Fj0rtoft argument [25j . we expect to have a direct cascade of energy with a 
£—17/5 sp ec trum flowing to large k and an inverse cascade of wave action with spectrum fc~ 3 flowing to small k. 



NUMERICAL RESULTS FOR TRUNCATED-LIA 

In the following we will consider numerical simulations of the system (|6j) under the conditions in which a stationary 
turbulent cascade develops. Energy and wave action are injected in the vortex filament by a white-in-timc external 
forcing (f)(z, t) acting on a narrow band of wavenumbers around a given kf. In order to have a stationary cascade, we 
need additional terms which remove H and N at small and large scales. The equation of motion obtained from ^ is 
therefore modified as 



w 



„'|2 



(~l) p vV 2p w - aw + 4>. 



(26) 



In (26) the small scale dissipative term (with p > 1) physically represents the radiation of phonons (at a rate 



proportional to v) and the large scale damping term can be interpreted as the friction induced by normal fluid at a 
rate a. 



Assuming the spectra ( 24 1 and ( 25 1 , a simple dimensional analysis gives the IR and UV cutoff induced by the 
dissipative terms. The direct cascade is removed at a scale k v ~ j/- 5 /( 10 p- 2 ) while the inverse cascade is stopped at 
k a ~ a 1 / 2 . Therefore, in an idealized realization of infinite resolution one would obtain a double cascade by keeping 
kf = 0(1) and letting v 1 a — > 0. In order to have an extended inertial range, in finite resolution numerical simulations 
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we will restrict ourselves to resolve a single cascade at a time by putting either kf ~ k a or kf ~ k v for the direct and 
inverse cascades respectively. 



We have developed a numerical code which integrates the equation of motion ( 26 ) by means of a pseudospectral 
method for a periodic vortex filament of length lix with a resolution of M points. The linear and dissipative terms 
are integrated explicitly while the nonlinear term is solved by a second-order Runge-Kutta time scheme. The vortex 
filament is initially a straight line (w(z,t = 0) = 0) and long time integration is performed until a stationary regime 
(indicated by the values of H and N) is reached. The ratio between the two terms in the series Q is H1/H2 — 20, 
confirming a posteriori the validity of the perturbative series ^ and the condition of the small amplitude perturbation 
e <C 1 in the derivation of equation ^ . 

The first set of simulations is devoted to the study of the direct cascade. Energy fluctuations are injected at a 
forcing wavenumber kf ~ 2 and the friction coefficient a is set in order to have k a ~ kf. Energy is removed at small 
scales by hyperviscosity of order p = 4 which restricts the range of dissipation to the wavenumber in a close vicinity 




FIG. 1: Wavenumber spectrum rik for a simulation of the direct cascade in stationary conditions at resolution M = 2048. 
Forcing is restricted to a range of wavenumbers 1 < kf < 3 and dissipation by phonon emission is modeled with hyperviscosity 
of order p = 4. The straight line represents the kinetic equation prediction ~ fc _17//5 . The inset shows the spectrum 
compensated with the theoretical prediction. 



In Figure [T] we plot the wave action spectrum for the direct cascade simulation, averaged over time in stationary 
conditions. A well developed power law spectrum very close to prediction (24 1 is observed over more than one decade 
(see inset). This spectrum confirms the existence of non-trivial dynamics with six- wave processes for the truncated 
Hamiltonian ([6]). 

The direct cascade of the full Biot-Savart Hamiltonian ^ was discussed by Kozik and Svistunov who gave the 
dimensional prediction ( 24 1 [1] , who later performed a numerical simulation of the nonlocal BSE to confirm the scaling 
0. 

We now turn to the simulation for the inverse cascade regime. To obtain an inverse cascade, forcing is concentrated 
at small scales, here kf = 300. In order to avoid finite size effects and accumulation at the largest scale [26], the friction 
coefficient is chosen in such a way that wave action is removed at a scale k a ~ 10. Figure [2] shows the spectrum for 
this inverse cascade in stationary conditions. In the compensated plot, a small deviation from the power-law scaling at 
small wavenumbers is observed, probably due to the presence of condensation ("inverse bottleneck") . Nevertheless, 
a clear scaling compatible with the dimensional analysis of the kinetic equation is observed over about a decade. 



DIFFERENTIAL APPROXIMATION MODEL 



Differential Approximation Models (DAMs) have proved to be a very useful tool in the analysis of fluid dynamical 
and wave turbulence [HI [16j [17l [T8j HH [20j [21] . These equations are constructed using a differential closure, such 
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FIG. 2: Wavenumber spectrum nk for a simulation of the inverse cascade in stationary conditions at resolution M — 2048. 
Forcing is restricted to a range of wavenumbers around kf = 300 and dissipation by phonon emission is modelled with 
hyperviscosity of order p — 4. The straight line represents the kinetic equation prediction ~ k~ 3 . The inset shows the 
spectrum compensated with the theoretical prediction. 



that the main scalings of the original closure (kinetic equation in our case) are preserved. In addition, there exist a 
family of simpler or reduced DAMs for which rigourous analysis can be performed upon their solutons. These appears 
to be quite helpful when the full details about the dynamics are not needed. Moreover, due to the DAMs simplicity, 
one can add physically relevant forcing and dissipative terms to the models. 

For the Kelvin wave spectra ( 24 1 and ( 25 1 , including the thermodynamical equipartition solutions ( 23 ) , the corre- 
sponding DAM is [5] 



~dt 



C_ 



,1/2 



duj 2 



n,V 1/2 



d 2 i 

duj 2 



(27) 



where k is the vortex line circulation, C is a dimensionless constant and u> = uj{k) = k/c 2 /47t is the Kelvin wave 
frequency. Notice, that the DAM is written in terms of frequency uj, rather than the wavenumber k, as in the case of 



the kinetic equation (21 1 



Equation ( 27 1 preserves the energy 



and wave action 



E 



N = 



-1/2, 



(28) 



(29) 



The forcing of Kelvin waves on quantized vortices arises from sharp cusps produced by vortex reconnections [S]. 
These reconnections can be interpreted in the DAM by the addition of the term [5] 



9n u 



dt 



(30) 



forcing 



The transfer of energy flux in the Kelvin wave turbulence proceeds towards high wavenumbers until, the Kelvin wave 
frequencies become large enough to excite phonons in the fluid and thus dissipate the energy of the Kelvin waves 
into the surrounding fluid. One can introduce sound dissipation derived from the theory of Lighthill in classical 
hydrodynamical turbulence [57] , and in the context of quantum turbulence has been done in [55] . This corresponds 
to the addition of the following term to the DAM [5] 



dn u 



dt 



= —vijjn, 



(31) 



radiation 
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(This expression is slightly corrected with respect to the relation ~ w 9 / 5 rt 2 of the one introduced dimensionally in [5 
to make it consistent with the more rigorous analysis of |28j). 

Finally, one can add an addition term that describes the effect of friction with the normal fluid component as 

follows, muni 



dt 



— —aujn u . (32) 

friction 



Thus, one may write a generalized DAM as 

= F(n^,uj) + \iv~ 2 - vuj^n 2 ~ aivn w . (33) 

where the nonlinear function F(n w , ui) is the interaction term, which in the case of the complete DAM is the RHS of 
equation ( p7| ). 

Other (reduced) DAMs include either, solutions for the direct and inverse cascade and no (thermo) equipartition 
solutions as in 

F( nuJ ,u) = LJ-V 2 d u (ntu 8 d u (u 3 / 2 n„)), (34) 

or just the direct energy cascade and the corresponding energy thermo solution (and no inverse cascade solutions): 

F(n u ,cj) = uj-^d^n^^d^Lon^)). (35) 

These have the advantage of only containing second order derivatives, and as such one may find analytical solutions 
for the steady state dynamics. For example, for model (34) we can ask the question, how does the vortex reconnection 
forcing build the energy flux in frequency space? For this, we leave the nonlinear transfer and the reconnection forcing 
terms and drop the dissipation term, and find a solution for the energy flux e w : 

e u = e - Acj- 1 / 2 (36) 

and the wave action density n^: 

= lOV^-V^^-i - hco-V 2 - Vo ) 1/5 (37) 

where eo and 770 are the asymptotic values of the energy and wave action fluxes respectively. 

For equation (35 1, without any forcing or dissipation terms, we find that the corresponding analytical steady state 
solution is 

n u ^Au-^^v + Blu- 1 , (38) 

which is a "warm cascade" solution, i.e. a direct energy cascade gradually transitioning into a thermalized "bottleneck" 
over a range of scales Qjj] . 



NUMERICAL RESULTS FOR DAMS 



We performed numerical simulations of the DAMs ( 27 1 , (34 1 and ( 35 1 , using a second order finite difference method. 
We set the resolution at M = 1024 points, while the phonon radiation dissipation term acts in the range u> > lo v — 800. 
The parameters A, v and the estimated asymptotic values of the fluxes eo, ?7o are listed in Table 1. The factor C/k 10 



of equation ( 27 1 has been fixed to unity. 



Initially, we simulate both the complete DAM (27 1 and the reduced DAM (34 1, with forcing, without friction, but 



with and without the dissipative term for phonon radiation. The results for the energy flux and the spectrum for 
the complete DAM are shown in Figure [3] (results for the reduced DAM are nearly identical and, therefore, are not 
shown) . The top panel in Figure [3] shows the energy flux. We see a good agreement with the analytical prediction 
( 36 1 over a large intermediate range of scales. In the case without the phonon dissipation, the numerical result for the 



energy flux follows perfectly the analytical prediction at high frequencies. In the case with the phonon dissipation, the 



agreement with the analytical prediction (36 1 is good in a long range up to very high frequencies, where the phonon 
dissipation suddenly kicks in. Such a sudden onset of dissipation is due to the abrupt growth of the phonon radiation 
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Reduced DAM 


Complete DAM 


Resolution 
Forcing amplitude 
Viscosity 
Energy flux 
Waveaction flux 


M = 1024 
A = 2 x 10" 9 
c = 2x 1(T 12 
e ^ 1.66 x 10" 9 
770 1.8 x 1(T 12 


M = 1024 
A = 1 x 10" 9 
u = lx 10~ 12 
e ~ 0.66 x 10~ 9 
770 ~ 0.75 x 10" 12 



TABLE I: Parameters and fluxes of the simulations for DAM models with v 7^ 



term as a function of the frequency. The wave action spectra are shown in the bottom panel of Figure [3] We see 
a good agreement with the analytical solution ( 37 1 where 770 is taken to be zero (flux of wave action 779 could only 
be generated by an extra forcing at the high-frequency end, which is absent in our case). We observe only a slightly 



steeper spectrum of n u 



-1.73 



compared with the analytical prediction of the direct energy cascade n u 



-17/10 



This agreement is remarkable because the solution (34) is strictly valid only for the reduced and not the complete 
DAM. This shows that the reduced DAM does pretty well in predicting the behaviour of the more complete nonlinear 
model. Naturally, in the case with the phonon dissipation, we sec a deviation from the analytical solution at very 
small scales (a rather sharp cut-off) . 




FIG. 3: Complete DAM: the energy flux e (top) and the spectrum (bottom), compared with the analytical predictions (36 1 and 
( |37[ l respectively. In each picture the results of the simulations with and without phonon dissipation are shown (the phonon 
dissipation can be seen as an abrupt cutoff of the flux and the spectrum). The spectrum is very close to the prediction of the 
reduced model (34 1: n u ~ u>~ 1,73 compared with ~ w - which is expected in the inertial range (10 < ui < 800). 
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In Figure [1] we show the effect of switching on the friction term (a ^ 0) for the complete DAM (again, results 
for the reduced DAM are very similar and are not shown) . We see that the presence of the friction has the effect of 
reducing the energy flux in a large region from frequencies of around wkIO upwards. Although the flux is reduced, 
we see that the spectrum has only a slight deviation from the predicted slope ( 37 1 . 

Finally we consider another DAM model, (35 1, which contains only the direct and the thermodynamical bottleneck 
(or warm cascade) of energy. We force the system as usual, however, the simulation is composed of two phases: the 
first is to get a direct cascade steady state, then lowering the viscosity v (i.e. dissipation at small scales) and then 
turning on the reflecting energy flux boundary condition at the smallest scale lu = uj max — M, after which the 
system evolves to a secondary steady state ( 38 1 . In Figure [5] we see a clear transition from the KZ solution towards 
the thermalized spectrum. The bottom panel is a zoomed in section of the crossover using a compensated spectrum, 
so that one can clearly make the distinction between the two power laws. We see a good agreement with the predicted 
power law behaviour of the KZ solution and of the thermalized solution. 




FIG. 4: Complete DAM with a/0 and v — or v 7^ 0. The energy flux is reduced at high w's, but the spectrum slope remains 
as in the non-dissipative case for a longer frequency range. 
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FIG. 5: Bottleneck effect for the second order model [5]. The spectrum (n(u) w.r.t. uj) compared with the two different power 
law predictions in the two ranges and the local slope analysis in the inset (a w.r.t. uS) which points out the transition from the 
rtuj ~ lo~ 1,7 behaviour to the n u ~ <^ _1 one. In the picture below, we show a zoom of the region where the bottleneck effect 
appears in a sort of compensated spectrum, nl^uj 27 ^ 2 on the y-axis and w on the x-axis. 



CONCLUSIONS 



In summary, we have introduced and studied various reduced models for Kelvin wave turbulence. Firstly, we 
introduced a truncated-LIA model for Kelvin wave turbulence, which we have shown to exhibit the same scalings 
and dynamical features present in the conventional BSE. We have used this model for numerical simulations of the 
direct and the inverse cascades and found spectra which are in very good agreement with the predictions of the 
WT theory. Secondly, we discussed differential approximation models, and introduced three such models to describe 
various settings of Kelvin wave turbulence, such as the direct energy cascade generated by vortex reconnections and 
dissipated via phonon radiation or/and the mutual friction with the normal liquid, as well as the bottleneck effect 
when the energy flux is reflected from the smallest scale. We performed numerical simulations of these cases, which 
showed good agreements with the predicted analytical solutions. 
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APPENDIX A - INTERACTION COEFFICIENTS IN BIOT-SAVART MODEL 

In this Appendix we review and extend the work of Kozik and Svistunov on the Kelvin wave cascade (KS04) 0] . 
They considered the full Biot-Savart Hamiltonian Q in 2D, and simplified the denominator by Taylor expansion. The 
criterion for Kelvin- wave turbulence is that the wave amplitude is small compared to wavelength, this is formulated 
as: 



£ (z 1 ,z 2 )= K f ) ' W ^ 2)l «l. 
\z\ - z 2 \ 



(39) 



KS04 find the Biot-Savart Hamiltonian ([2| expanded in powers of e (H = Ho + H\ + H 2 + H3, here Ho is just a 
number and is ignored) is represented as: 



Hi 
H 2 
H 3 



K f dz\dzi 
8tt J \zi - z 2 \ 
k f dz\dz 2 
327T J \zi - z 2 
k f dz\dz 2 



2Re 



3e 4 - Ae 2 Re 



647T J \zi~ z 2 \ 



6e 4 i?e 



(w *(zi)w (z 2 ) 

(w *{zi)w (z 2 ) 
*(zi)w (z 2 ) 



(40) 



One would like to deal with the wave-interaction Hamiltonian (40 1 in Fourier space by introducing the wave 
amplitude variable Wk- Using the Fourier representation for variables w(zi) — J dkwke lkzi and w(z 2 ) — J dkwke tkZ2 in 
equations (40), one introduces more integration variables, the wavenumbers. Moreover, invoking a cutoff at a < \z\—z 2 \ 
because of the singularity present in the Biot-Savart Hamiltonian (p| as \zi — z%\ — > 0, KS04 derived the coefficients 



of Hi, H 2 and H3 in terms of cosines in Fourier space [4]. Once the equations (40 1 are written in terms of wave 
amplitudes, one can define variables z_ and z + as variables z_ = \zi — z%\ which ranges from a to 00 and z + = z\ + z 2 
ranging from —00 to 00. One can then decomposes all variables of type Zi and z 2 into variables z + and z_. The 
cosine functions arise due to the collection of exponentials with powers in variable z_ . Subsequently, the remaining 
exponentials with powers of variable z + can be integrated out w.r.t. z + , yielding the corresponding delta function for 
the conservation of wavenumbers. The explicit formulae of the four-wave W1234, and the six-wave T123456 interaction 
coefficients derived by Kozik and Svistunov can be written as: 



W1234 

T'l23456 



2tt 

K 

16tt 

K 

16tt 



[A-B], 



[§D - E] , 
[3P - 5Q] 



(41) 
(42) 
(43) 



where A, B, D, E, P and Q are integrals of cosines: 
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.4 
B 
D 
E 

P 

Q 



dz- 



■k 2 C k , 



00 



[I-C1-C2- C 3 



C 4 



cl + c 43 + c 2 4 l 



[feifc 4 (C 4 + Ci - C 43 - C 2 4 ) + kika (C 3 + Ci - C 



43 



c 3 



fc 3 fc 2 (c 

T5~ 



c 2 -c 



43 



Cf) + k 4 k 2 (C 4 + C 2 - C 43 - C 3 )], 



fcefc2[C 2 — Cf — C23 + cf 3 — C 2 + C 2 5 + C 23 — + c 6 



c 



56 



-C3 6 



cf - c 



46 



C 456 + Cf -C12], 

P°° d? 

\ ^[i~c 4 -c 1 + c 4 -c 6 + c 46 + cf - c 46 - c 5 

+C 65 - C 456 - Cf 6 + C; 



c 



45 



Ci b - cf 



45 



'23 



C3 + C3 + C13 — C 13 + Cg 



c 



46 



13 "T" °2 



+c| - c 3 45 



C5 1 s~i6 /~y56 1 s~i 1 /t4 /" ( 
13 "r ^2 — ^3 ~r ^12 t U 2 ~ U2 



(44) 
(45) 
(46) 

(47) 
(48) 

(49) 



where the variable, z_ = j^i — z%\ and the expressions C, are cosine functions such that Ci = cos(fciz_), C\ — 
cos((fc 4 — ki)z-), Cf 5 = cos((k 4 + feg — C 4 | = cos((/c4 + k 5 — k\ — k 2 )z-) and so on. 



We integrate the Fourier representations of Hi, H 2 and H 3 , namely equations (44 1, (451, (46 1, (|47|, (48 1 and (49 1 
using integration by parts, and apply the following cosine identity |29j . 



cos(z) 



dz 



r cos(z) - i 

-7 — ln(a) — / dz 



= -7 - ln(a) - ^ 



j^2fc(2fe)! 
- 7 -ln(a) + C(a 2 ). 



(50) 



Neglecting terms of order ~ a and higher we calculate the frequency and interaction coefficients of equations ( 40 1 



Hi 
H 2 



LUk\wk\ 2 dk, 

dk 1234 W 1234 wl w 2 w 3 w 4 SH, 



H 3 = I rffci23456Ci23456^*^2 u; 3 u; 4^5'W6^456) 



where 



L0 k 



w 1234 



-k 2 



C 



4tt 

K 

64?r 

K 



hi 



k cB a 



Kk 2 1 / k 

47T \ fc eff 



k x k 2 k 3 k 4 



1 + 47 - 4 In 



1 



ke&a 



-^1234) 



123456 



1287T 



k 1 k 2 k 3 k 4 k 5 k 6 



1 - 47 + 4 In 



fc c ff a 



G 



123456- 



(51) 
(52) 
(53) 

(54) 
(55) 
(56) 



We use the notation that fc e ff is the mean value of wavenumbers, 7 = 0.5772 ... is the Euler constant and -F1234 and 
G123456 are logarithmic terms of order one that are shown below, 



F 



1234 — 



iV 4 , f N 

, In 

16tt L ^ 24 V fccff 
Nei v 



NeJ 



fccff 



(57) 



Gi 



23456 



EN 4 
kfiky In 
6 2 24 

NeK 



L6tt L 



N 

fccff 



720 



N 



16 

(58) 



7 
J 



L = 



{-[i}-[2],-[%-[%[%nti}} 

{{[ 4 ],[i],-[ 43 ],-[!]} i=4J=1 , 

{pi.w.-n.-B]}^. 

{[ 4 ],[2],-[ 43 ],-[ 4 ]} i=4J=2 }, 

{[ 2 ],-g],-[23],g 3 ]r 



[iMfLbU-IiU 6 ] 



r5(il 



[i 6 ],-n,n,[f],- 



{ 



i}. 



[Mi],- 



r456i 



[13]) [3]) - [3 5 L - [l3]) S) _ [ 3 5 ]) [12], [2]) ■ 



[13J 
r65i 



[ 5 ],[ 45 ],[?],-[f], 



[13]' [3]' b 6 ]' 



(59) 



(60) 



(61) 



]}• (62) 
The notation used for the logarithmic terms is as follows, for TV e K — =b[£ g ], the corresponding term in G123456 is 



3k, (kr) H~ kn kr) ( k< 



kg k r 



(63) 



We checked numerically that after the canonical transformation of the four-wave dynamics, that in the full six-wave 
interaction coefficient C123456, the LIA contribution, which is of order A(fc e ff) = hi(fc c ffa) in the interaction coefficient 
(42 1 directly cancels with the LIA contribution from the six- wave interaction in d43|. The canonical transformation, 



which is explicitly written in the main text, equation ( 14 1 can be expanded by the means of a small parameter 1 / K(k c s ) . 



We expand the denominator, which is just a function of the linear frequency around the A(k c g) contribution, 



0J k - k 2 [A(fceff) + A] , 

1 1 A 

LO k ~ A(fc cff ) A(fceff) 2 

Using this strategy, we can represent the final six-wave contribution in terms of powers of 1/A(fc c ff). We find in this 
case that the leading term of order A(fc c ff), which corresponds to the local dynamics drops out, and the next leading 
order is of order unity. Moreover, we find that this leading contribution to the six-wave interaction coefficient, which 
is A(fc e fj) independent is also k e s independent. 

Unfortunately, there is no magic cancelations or simplifications in the resulting six-wave coefficient and it remains 
extremely large and complicated even after the integrations and the asymptotical limits of large A(k c g) taken in this 
Appendix. This is why we think that the truncated-LIA model, where the six-wave coefficient is simply the product 
of the six wavenumbcrs, is helpful for understanding the basic properties of Kelvin wave turbulence, before a more 
complete model based on the BSE can be attacked. However, we think that the derivations made in this Appendix 
will help us to at least establish whether or not the KZ spectra of Kelvin wave turbulence are local or nonlocal, and 
a work is underway in this direction. 



(64) 
(65) 



APPENDIX B - INTEGRABILITY OF 2D-LIA 

We show the details of the derivation that the LIA of the BSE yields the 2D-LIA model Q, i.e. that the cutoff 
operation commutes with making the 2D reduction. In view of the integrability of the original LIA, this amounts to 
a proof of integrability of the 2D-LIA model Q . 

The LIA of the BSE can be written as 
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r = /3r' x r", 



(66) 



1/2 

where the notation for the differentiation operator ' is d/dl, and where I = (l + |w| 2 ) is the arc length. The 
two-dimensional representation of a vortex line can be described by a vector r = zz + w. The vector w = (x(z), y{z)) 
being a function of z, orientated in the xy-plane. 

Applying the chain rule to rewrite all derivatives to be with respect to z (from this point on ' will refer to d/dz), 
LIA can be written as 



r = (3 



1 



|w'| 2 r 3/2 (h w" + w' xw") 



(67) 



With a little geometrical intuition [2], one can show w = r — (r ■ z)(z + w'). In addition, both w' and w" are 
perpendicular to z direction, thus, one can represent w' x w" = ((w' x w") ■ z)z = Az. Then equation (67) can be 
reduced to 



W = /3(l + |w'| 2 ) 3/2 [z x w" - AW] 



We will now show that equation Q is equivalent to equation (68 1. First, we must change our representation of 
from a complex variable w{z) = x(z) + iy(z) to vector notation w = (x(z), y(z)). Equation Q is equivalent to 



(68) 
w 



w 



1. d [ 
w = t: z x — — . z 
2 dz I ^/l + | w f 



(69) 



Expanding, keeping track of (3 and applying the vector identity (a x b) x c = (c • a)b — (c • b)a, equation (69 1 can be 
rewritten as 



w = /3(l + |w'| 2 ) 3/2 [z x w" - AW] 



(70) 



This is exactly the same result as Equation (68). The 2D-LIA model Q is equivalent to the LIA of the BSE (66 1 
and so the 2D-LIA model is indeed integrable. We remind, however, that both models Q and (66 1 can only be 
used continuously until the moment when they predict self-crossings of vortex lines (absent in weak turbulence), and 
therefore the integrability of these models may only be used for predicting the vortex line motion in between of the 
reconnection events (see the footnote on page 4). 
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